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Abstract 

To retrieve temperature and humidity profiles from SSM/T and AMSU, it is important to 
quantify the contribution of the Earth surface emission. So far, no global estimates of the land 
surface emissivities are available at SSM/T and AMSU frequencies and scanning conditions. 
The land surface emissivities have been previously calculated for the globe from the SSM/1 
conical scanner between 19 and 85 GHz. To analyze the feasibility of deriving SSM/T and 
AMSU land surface emissivities from SSM/I emissivities, the spectral and angular variations of 
the emissivities are studied, with the help of ground-based measurements, models and satellite 
estimates. Up to 100 GHz, for snow and ice free areas, the SSM/T and AMSU emissivities can 
be derived with useful accuracy from the SSM/I emissivities. The emissivities can be linearly 
interpolated in frequency. Based on ground-based emissivity measurements of various surface 
types, a simple model is proposed to estimate SSM/T and AMSU emissivities for all zenith angles 
knowing only the emissivities for the vertical and horizontal polarizations at 53° zenith angle. 
The method is tested on the SSM/T-2 91.655 GHz channels. The mean difference between the 
SSM/T-2 and SSM/l-derived emissivities is < 0.01 for all zenith angles with an r.m.s. difference 
of % 0.02. Above 100 GHz, preliminary results are presented at 150 GHz, based on SSM/T-2 
observations and are compared with the very few’ estimations available in the literature. 

Keywords 

Microwave radiometry, surface emissivity, ATOVS. 


I. Introduction 

The Special Sensor Microwave/Temperature 1 and 2 (SSM/T-1 and -2) and the Ad- 
vanced Microwave Sounding Units A and B (AMSU-A and -B) are both cross-track tem- 
perature and water vapor profilers with similar frequencies, but AMSU has better spatial 
resolution. The SSM/T instruments are on board the Defence Meteorological Satellite 
Program (DMSP) polar orbiting satellites. SSM/T-1 has 7 channels in the 0 2 absorption 
band around 60 GHz for temperature sounding of the atmosphere [1]. The SSM/T-2 is a 
water vapor profiler with 5 channels, three in the H 2 0 absorption line at 183.3 GHz and 
two window channels at 91.655 and 150 GHz [2]. The Advanced Microwave Sounding Unit 
(AMSU), part of the Advanced TIROS Operational Vertical Sounder (ATOVS), replaces 
the Microwave Sounding Unit (MSU) on the previous NOAA polar orbiters. AMSU in- 
cludes a temperature sounder (AMSU-A) with 15 channels, most of them located in the 
0 2 absorption band around 60 GHz, and a humidity sounder (AMSU-B) w T ith channels 
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Microwave land surface emissivities over the globe have been estimated from the Special 
Sensor Microwave/ Imager (SSM/I) observations, by removing the contributions of the at- 
mosphere, clouds and rain using ancillary satellite data [16]-[17]. The SSM/I instrument 
is described in [18]. Cloud-free SSM/I observations are first isolated with the help of col- 
located visible/infrared satellite observations (International Satellite Cloud Climatology 
Project (ISCCP) data [19], Then, the cloud-free atmospheric contribution is calculated 
from an estimate of the local atmospheric temperature-humidity profile (National Centers 
for Environmental Prediction (NCEP) analysis [20]. Finally, with a surface skin temper- 
ature derived from IR observations (ISCCP estimate), the surface emissivity is calculated 
for all the SSM/I channels. The emissivities are estimated for a 53° observation angle at 
19.35. 22.235, 37.0 and 85.5 GHz for both vertical and horizontal polarizations with the 
exception of 22 GHz, which has vertical polarization only. The emissivities are available 
on a 1/4° grid, compatible with the ~ 30 x 30km ISCCP DX grid and with the SSM/I 
observations which are sampled at 25 km. The standard deviation of the day-to-day varia- 
tions of the retrieved emissivities within a month is typically about 0.012 for all the SSM/I 
frequencies, which is an upper limit on uncertainty of these estimates. Biases arising from 
uncertainties in the IR emissivity are <0.02. 

Similar technics could be applied to SSM/T and AMSU observations to derive the land 
surface emissivities for each frequency and scanning angle, but because the viewing angles 
of SSM/T and AMSU are not constant, most scenes on the globe are not seen more than 
once a month under clear sky conditions with the same angle. Thus, to obtain an adequate 
climatology, a long time series of data (3 years at least) would have to be processed before 
having reliable estimates of the natural variability of the surface emissivities. AMSU 
is now operational since January 1999, and as the data become available, calculations 
of the emissivities will be performed at the Centre de Meteorology Spatial at Lannion 
(France) with a method similar to the one developed for SSM/I. However, before an 
adequate time series of the emissivities become available, some practical alternative has 
to be implemented now in order to efficiently process the satellite data over land. 

We examine the feasibility of estimating the SSM/T and AMSU emissivities over the 
globe from the previously retrieved SSM/I emissivities, taking into account the different 
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observations. The p — polarized soil emissivity ts p is given by: 

ts, = i - in + 47(<?s - 1) - “•''*(»>) a) 

where p and q stand for the vertical or horizontal polarizations, t S s' c is the p — polarized 
specular emissivity; h$ and Qs are respectively the roughness and the polarization mixing 
parameters; N$ is an exponent fitted to reproduce the angular variations. The soil dielec- 
tric permittivity is calculated from the model described in Calvet et al. [22]. The values 
of the average best fit parameters hs, Qs and Ns are given in Table 2 for the three plots 
and the three frequencies. 

To account for the effect of the vegetation canopies, the emissivity cc of a vegetated 
surface is computed as follows, for the p— polarization: 

£cp = ^Sp x fsp + (1 fs P ) x fv (2) 

where es P is the emissivity of the bare soil calculated from the previous equations, ty is 
the emissivity of the vegetation cover, and fs p is the fraction of bare soil. 

The modeling is based on several assumptions. First, it is assumed that the vegetation 
emissivity simply depends on the crop type and frequency and that it does not depend 
on polarization and incidence angle. Second, since crops are generally arranged in row’s, 
geometrical effects may be significant. To account for this effect a model is developed for 
polarization p: 

fsp = (1 -QpX tan(0)) x f so (3) 

where q p is a best fit coefficient and fso is the fractional coverage of soil seen at nadir. 
However, this last equation can be considered as a modeling refinement for crop covers 
w’ith row structure and is likely to be unnecessary for large scale footprints which include 
a variety of surface types. 

The simple model described here for vegetation and soil provides a good fit to the avail- 
able observations w’ith an average r.m.s. error of 0.002 between measured and simulated 
emissivities. This model is used in our study to analyze the frequency and angular depen- 
dence of the land surface emissivities up to 100 GHz. 
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The RADTRAN surface mode] calculates vertically and horizontally polarized surface 
emissivities for various surface types for frequencies up to 40 GHz [24], The modeling 
approach for vegetation and bare soil is based on radiative transfer theory where the 
vegetation is treated as layers of continuous random media bounded by an underlying 
homogeneous soil layer. The RADTRAN model predicts an increase in emissivity with 
frequency for vegetation and for bare soil. However, the comparison between simulations 
and SSM/I data provided in [24] shows that the surface emissivity is somewhat overesti- 
mated by the model at 37 GHz (case c in Figure 3 in the paper). 

B. From satellite-based estimates 

The SSM/I frequencies range from 19 to 85 GHz. For most surfaces, the emissivities in 
this frequency range vary smoothly with frequency for both orthogonal polarizations at 
53 u incidence angle [16]. For nine vegetation classes derived from Matthews [25], Figure 2 
shows the frequency dependence of the mean surface emissivity at 53° zenith angle for each 
surface type, as calculated for the Meteosat area (Africa plus large portions of Europe and 
Western Asia) for October 1991. Whatever the vegetation type, the emissivities slowly 
decrease with frequency for both orthogonal polarizations. For the vertical polarization, 
the emissivity change between 19 and 85 GHz rarely exceeds 0.05 and is smaller over dense 
vegetation than over bare soil; for the horizontal polarization, the changes are smaller still 
(always < 0.025). Under snow and ice conditions, the surface emissivity varies more 
quickly with frequency and these surface types will have to be studied further. 

Alternative estimates of microwave surface emissivities from satellite in this frequency- 
range are scarce. However, all the available estimates from SSM/I observations show that, 
for various surfaces (bare soil, vegetated soil), the surface emissivity decreases with increas- 
ing frequency for both vertical and horizontal polarizations at 53° incidence. Choudhury 
[26] analyzed rain forest and desert locations for the period January 1988 to December 
1989 and found that emissivities at 37 GHz are almost always lower than at 19 GHz for 
both orthogonal polarizations. For a 70-day period over the centra] United States, Jones 
and Vonder Haar [27] also observed a decrease in emissivity with increasing frequencies 
from SSM/I measurements at both polarizations. The same trend is confirmed by Xiang 
and Smith [28] with SSM/1 observations of the Sahelian region. 
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IV. The angular dependence of the microwave land surface emissivities 

The SSM/T and AMSU instruments are both cross-track scanners. SSM/T has 7 scan 
positions 9 S from -39° to +39° and SSM/T-2 has 28 positions from -40.5° to +40.5°. These 
are satellite view angles which translate into local zenith angles 9 Z up to 47.4° near the 
edge of the scan due to the curvature of the Earth. AMSU- A has 30 scan positions at 3.3° 
intervals from -14.5x3.3° to -1-14.5x3. 3° while AMSU-B has 90 positions at 1.1° intervals 
from -44.5x1.1° to -1-44.5x1.1° which translate into local zenith angles 9 Z up to 58.5°. The 
polarization measured by SSM/T and AMSU rotates with scan angle due to the rotating- 
reflector/fixed-feed type of antenna design. If 6 S is the scan angle and 0 Z is the local zenith 
angle, then the SSM/T or AMSU surface emissivity t(9 z ) seen for a local zenith angle 9. 
is given by; 

e{9 z ) = t p {9,) cos 2 (0 s ) -|- £,(0 z )sin 2 (0 s ) (4) 

c p{9 : ) and t q {9 : ) are the two orthogonal polarized surface emissivities at 9 Z local zenith 
angle. Depending on the channels, p will represent the vertical or the horizontal polariza- 
tion. The polarization p seen when the incidence is close to nadir (i. e. for 9 Z =9 $ very 
close to 0°) is indicated for each channel on Table 1. 0 4 =45° corresponds to 9 Z = 53° 
(which is also the SSM/I zenith angle); for this angle, c(53°) = (c p (53°) + e,(53°))/2. 

The polarization state for SSM/T-2 is sometimes given as “unspecified”. In some studies 
([29] for instance). SSM/T-2 has been assumed to observe vertical polarization at nadir 
whereas horizontal polarization at nadir has been assumed by Wessel and Boucher [30] 
in their comparison of the SSM/I and SSM/T-2 window channels near 90 GHz. From 
comparisons between observations and simulations, Burns et al. [31] concluded that the 
instrument is observing the horizontal polarization at nadir and this has been confirmed 
to them by information from the Aerojet system engineers for the SSM/T-2 project. 

The SSM/I emissivities are only available at one zenith angle (53°) and for two orthog- 
onal polarizations, giving no information on the emissivity angular dependence. 

Model results derived from the INRA measurements are analyzed to estimate the angular 
dependence of the emissivity. The INRA model is used to simulate the surface emissivi- 
ties at frequencies up to 100 GHz, for angles between 0° and 60° and for two orthogonal 
polarizations. The results are presented in Figure 3 for different surface types (smooth 
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to the cross-track scanners, a single set of coefficients a„ is calculated to minimize the r.m.s. 
errors between the fitted F(ey(53°), e#(53°), 9 Z ) and the simulated emissivities (solid lines 
in Figure 3) for all the cases presented in Figure 3. The fit at 90 GHz is added in Figure 3 
for each surface type, along with the r.m.s. error between the fit and the model. The r.m.s. 
error is <0.015, whatever the surface type. The ’x’ symbol in Figure 3 indicates the mean 
of the two orthogonal polarizations at 53°: The r.m.s. error resulting from the use of this 
mean value regardless of the zenith angle is indicated in brackets. The r.m.s. error given 
by the angular dependent function is better or comparable to the r.m.s. error obtained 
with a fixed mean value. The corresponding a n coefficients are indicated on Table 3 for 
23.8. 36.5 and 90 GHz. for the two possible polarization patterns, i. e. with horizontal or 
vertical polarizations close to nadir. For each a n coefficient, a linear regression in frequency 
is applied and the corresponding a n (f) are also presented in Table 3 as a function of the 
frequency /. This is equivalent to applying a linear regression to the whole function. 

The function in equation 5 is then tested on different surface types using the RADTRAN 
model and the a n (f) previously calculated. The results are presented in Figure 4 for two 
AMSU-A frequencies. For bare soils and for vegetated areas, the function in equation 
5 represents well the model angular dependence, with r.m.s. errors lower than 0.015. 
For open water surfaces, the angular dependence of the surface emissivity is not well 
reproduced. These surfaces are characterized by low emissivities at 53° for the horizontal 
polarization (<0.5 at 19 GHz): They represent less than 1% of the surface emissivities as 
calculated from SSM/I and are concentrated in lake or coast areas. These pixels can be 
processed separately. 


V. Verification of the method for SSM/T-2 observations at 91.655 GHz 

AND AT FREQUENCIES ABOVE 100 GHZ 

Using a linear interpolation in frequency and the angle-dependent model described above, 
SSM/T and AMSU emissivities up to 100 GHz can be estimated from SSM/I emissivities. 
We first verify this method against SSM/T observations at 91.655 GHz and then we 
examine a possible extension of the method above 100 GHz. 
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SSM/T-2 emissivities at 91.655 GHz and the SSM/I derived emissivities are presented 
on Figure 5 for three zenith angles. The results obtained when using a mean emissivity 
calculated at 53°, regardless of the zenith angle, are also shown with dashed lines. When 
using the angle-dependent model, the mean error (<0.01) and the r.m.s (ss0.02) are similar 
or lower for all scanning angles than when using the mean emissivity. Note that more than 
half the r.m.s. differences in Figure 5 can be accounted for by the intrinsic variations in the 
SSM/I emissivity values; this variability appears in a single month and can be expected 
to appear year-to-year [16]. 

B. Possible extension of the method to frequencies above 100 GHz ? 

Figure 6 presents the 150 GHz emissivities for February 1995 for three zenith angles, 
along with the emissivity difference between 150 and 91.655 GHz. For desert and sparse 
vegetation, the emissivities at 150 GHz are lower than at 91.655 GHz. This is not always 
the case for densely vegetated areas. These results are consistent with Felde and Pickle 
[29]. who also found different spectral variations depending on the surface type. However, 
as previously mentioned, Hewison [15] found higher emissivity at 150 GHz than at 89 GHz 
at nadir for all the surfaces they observed. Several 5° x5° areas, where the surface types are 
supposed to be homogeneous, have been selected to illustrate the angular dependence of 
the surface emissivities at 91.655 and 150 GHz, as estimated from SSM/T-2 measurements 
(Figure 7). For each zenith angle, the mean value is indicated along with the associated 
standard deviation. As expected, at 91.655 GHz, the angular dependence for bare soil 
and sparse vegetation is stronger than for densely vegetated areas. Despite a large scatter 
in the results, it appears that the angular dependence of the emissivities over tropical 
forest is stronger at 150 GHz than at 91 GHz. Over deciduous forests, Hewison [15] also 
measured a larger than expected angular dependence at 150 GHz (larger than the angular 
dependence they observe at 89 GHz). At 150 GHz, the standard deviations of the results 
are larger especially in tropical areas for large zenith angles. This is probably related to 
the larger atmospheric contribution at this frequency, especially in tropical areas. 

Figure 8 shows the change in the estimated emissivity for a given change in the atmo- 
spheric contribution (increase in the water vapor continuum absorption or increase in the 
water vapor column abundance) at 91 and 150 GHz for three standard atmospheres. The 
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the emissiv ities can be linearly interpolated in frequency. The scanning and polarization 
patterns of SSM /T and AMSU are such that the angular dependence of the emissivity seen 
by these instruments is rather small. Based on ground-based emissivity measurements of 
various surface types, a simple model is proposed to estimate SSM/T and AMSU emis- 
sivities for all zenith angles knowing only the emissivities for the vertical and horizontal 
polarizations at 53° zenith angle. The method is tested on the SSM/T-2 91.655 GHz 
channel. The spatial variability of the surface emissivities is well captured. The mean 
difference between the SSM/T-2 emissivities and the SSM/I derived emissivities at 91.655 
GHz is lower than 0.01 for all zenith angles with an r.m.s. difference of ~ 0.02. An atlas 
of the SSM/I land surface emissivities is available with a 1/4° resolution. The inter- and 
intra-annual variability of the emissivities is now under study. With the method devel- 
oped in this paper, emissivity maps at AMSU-A frequencies and scanning conditions are 
being prepared for the French Meteorological Office in Lannion (Centre de Meteorologie 
Spatiale). They will be used as emissivity first guesses in the temperature profile retrieval 
scheme. 

Above 100 GHz, preliminary results have been presented at 150 GHz, based on SSM/T-2 
observations. These results are compared with the very few estimations available in the 
literature. It appears that a simple frequency extrapolation from 90 GHz to 150 GHz will 
not give satisfactory results. Longer time series of satellite data at 150 GHz will have to 
be processed, and because of the sensitivity of this frequency to water vapor absorption, 
special attention will have to be paid to both the water vapor profile estimation and the 
water vapor absorption model. 
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Table 1. SSM/T and AMSU characteristics 


lo 


Channel 

number 

Frequency 
( GHz ) * 

Polarization 
at nadir 

Resolution 
at nadir (km) 

Atmospheric 

transmission 

(tropical) 

Atmospheric 
transmission 
(winter subarctic) 

1 

50.500 

H 

SSM/T- 1 

175 

.62 

.67 

2 

53.200 

H 

175 

.20 

.22 

3 

54.350 

H 

175 

.02 

.02 

4 

54.900 

H 

175 

.00 

.00 

5 

58.400 

V 

175 

.00 

.00 

6 

58.825 

V 

175 

.00 

.00 

7 

59.400 

V 

175 

.00 

.00 

8 

91.655 

H 

SSM/T- 2 

88 

.60 

.91 

9 

150.000 

H 

54 

.23 

.84 

10 

183 . 31 ± 7.00 

H 

48 

.00 

.40 

11 

183 . 31 ± 3.00 

H 

48 

.00 

.07 

12 

1 83.31 ± 1 .00 

H 

48 

.00 

.00 

1 

23.8 

V 

AMSU- A 

50 

.78 

.99 

2 

31.4 

V 

50 

.89 

.96 

3 

50.3 

V 

50 

.63 

.68 

4 

52.8 

V 

50 

.29 

.32 

5 

53 . 596 ±.l 15 

H 

50 

.11 

.13 

6 

54.40 

H 

50 

.02 

.02 

7 

54.94 

V 

50 

.00 

.00 

8 

55.50 

H 

50 

.00 

.00 

9 

57 . 290 = i / 

H 

50 

.00 

.00 

10 

i /±.23 7 

H 

50 

.00 

.00 

11 

v ± .322 ± .048 

H 

50 

.00 

.00 

12 

v ± .322 ± .022 

H 

50 

.00 

.00 

13 

v ± .322 ± .010 

H 

50 

.00 

.00 

14 

u± .322 ± .0045 

H 

50 

.00 

.00 

15 

89.0 

V 

50 

.61 

.91 

16 

89 . 0±.9 

V 

AMSU-B 

15 

.61 

.91 

17 

150 . 0 ±. 9 

V 

15 

.23 

.84 

18 

183 . 31 ± 1.00 

V 

15 

.00 

.00 

19 

183 . 31 ± 3.00 

V 

15 

.00 

.07 

20 

183 . 31 ± 7.00 

V 

15 

.00 

.40 



51 


Table 2. Soil parameters for the INRA model 


Soil type 

Frequency (GHz) 

h SOIL 

Qsoil 

N SOIL 

Smooth sod 

23.8 

0.43 

0.23 

0.0 


36.5 

0.15 

0.33 

0.0 


90.0 

0.41 

0.31 

0.0 

Rough sod 

23.8 

0.57 

0.43 

0.5 


36.5 

0.62 

0.45 

0.5 


90.0 

0.71 

0.40 

0.5 

Wheat field 

23.8 

0.40 

0.37 

0.0 


36.5 

1.00 

0.37 

0.0 


90.0 

0.50 

0.40 

0.0 



ev EG 

Qv 

Qh 


23.8 

.965 

0.43 

0.23 


36.5 

.980 

0.08 

1.05 


90.0 

.980 

0.23 

0.23 



Table 3. The angular model parameters 


Frequency (GHz) 

ao 

ai 

a 2 

a 3 

23.8 

0.13 

V’ nadir 
-5.99e-3 

5.21e-4 

-0.86e-5 

36.5 

0.24 

-7.56e-3 

6.33e-4 

-].09e-5 

90.0 

0.37 

-9.46e-3 

7.61e-4 

-1.35e-5 

f 

3.27e-3x f 

-4.74e-5x f 

3.26e-6x f 

-0.66e-7 x f 


+0.08 

-5.29e-3 

+4.75e-4 

-0.77e-5 

23.8 
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